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Abstract 



Closed analytical expressions for scattering intensity and other global struc- 
ture factors are derived for a new solvable model of polydisperse sticky hard 
spheres. The starting point is the exact solution of the "mean spherical ap- 
proximation" for hard core plus Yukawa potentials, in the limit of infinite 
amplitude and vanishing range of the attractive tail, with their product re- 
maining constant. The choice of factorizable coupling (stickiness) parame- 
ters in the Yukawa term yields a simpler "dyadic structure" in the Fourier 
transform of the Baxter factor correlation function qij(r), with a remarkable 
simplification in all structure functions with respect to previous works. The 
effect of size and stickiness polydispersity is analyzed and numerical results 
are presented for two particular versions of the model: i) when all polydisperse 
particles have a single, size-independent, stickiness parameter, and ii) when 
the stickiness parameters are proportional to the diameters. The existence 
of two different regimes for the average structure factor, respectively above 
and below a generalized Boyle temperature which depends on size polydis- 
persity, is recognized and discussed. Because of its analycity and simplicity, 
the model may be useful in the interpretation of small-angle scattering exper- 
imental data for polydisperse colloidal fluids of neutral particles with surface 
adhesion. 
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I. INTRODUCTION 



A theoretical determination of scattering intensity and structure factors for fluids with a 
large number p of components is, in general, a difficult task. In particular, this is true when 
some kind of polydispersity is present, as occurs in colloidal or micellar solutionsEHl. Size 
polidispersity means that macroparticles of a same chemical species exhibit several different 
dimensions within a discrete or continuous set of possible values. Interaction polydisper- 
sity then denotes a similar, and usually correlated, dispersion of parameters (charges, etc.) 
defining the strength of interaction potentials. Even when all macroparticles belong to a 
unique chemical species, a polydisperse fluid must therefore be treated as a multi-component 
mixture, with very large p values (of order 10 1 4- 10 3 ) or, in the infinite-component limit 
p — ► oo, with an idealized continuous distribution of some properties^. 

In structural studies on polydisperse colloidal fluids, a key role can be played by the 
models for which the Ornstein-Zernike (OZ) integral equations of the liquid state! admit 
analytical solutions leading to closed expressions of scattering functions for any finite p and 
even for p —* oo. A sufficient condition for this is that the Fourier transforms %j (k) of the 
functions %(r), solutions of the Baxter factorized version of the OZ equations^, have a 
peculiar mathematical form, which we refer to as dyadic structure! and will be illustrated 
in Section II C. Using the dyadicity, the explicit inversion of a related p x p matrix Q (k) is 
always possible for arbitrary p and closed analytical expressions for the "partial" structure 
factors Sij(k^ can be obtained. The scattering intensity and other "global" structure factors 
are then calculated as weighted sums of all partial structure factors. 

Usually, the above sums are performed numerically by evaluating p(p + 1) /2 indepen- 
dent contributions Sij(k) at each /clliHll. This procedure becomes numerically demanding 
for large p, as required in polydisperse mixtures. On the other hand, we stress that the 
dyadicity property also enables an alternative route (followed in the present work) which 
avoids the explicit computation of individual Sij(k). The weighted sums can, in fact, be 
worked out analytically, by a procedure originally proposed by VrijEl and referred to as 
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"Vrij's summation" hereafter. The resulting closed analytical expressions of "global" scatter- 
ing functions hold true for any number p of components, can be easily applied to polydisperse 
fluids even in the limit p —>■ oo, and are particularly suitable to fit experimental scattering 
data. 

VrijH 

first obtained a closed expression for the scattering intensity of polydisperse 
hard spheres (HS) within the Percus-Yevick (PY) approximation. Gazzillo et a/.i derived 
similar formulas for polydisperse charged hard spheres (CHS), by using the corresponding 
analytical solution within the "mean spherical approximation" (MSA). 

In the present paper we extend the approach previously exploited for polydisperse HsEH 
and CHsi to polydisperse "sticky" hard spheres (SHS). This simple model adds to an inter- 
particle hard core repulsion an infinitely strong attraction at contact, and can be applied to 
real colloidal fluids of neutral spherical particles with a van der Waals (or dispersion) force 
of attraction, working at very short distances. 

BaxteiS proposed the one-component original version of this model, solved the OZ equa- 
tion with the PY closure and found that such a system presents a liquid-gas phase transition. 
The adhesive contribution in Baxter's Hamiltonian is defined by a particular limiting case 
of a square-well tail in which the depth goes to infinity as the width goes to zero, in such a 
way that the contribution to the second virial coefficient remains finite but not zero (Baxter 
limit). StellB found that SHS of equal diameter in the Baxter limit, when treated exactly 
rather than in the PY approximation, are not thermodynamically stable. Nevertheless, the 
PY solution for SHS as a useful colloidal model has received a continuously growing interest 
in the last two decades0@, partially because of its capability to exhibit a gas-liquid phase 
transition. 

Perram and Smith0 and Barboy and TenneS extended Baxter's work to p-component 
fluids, using the same kind of Hamiltonian and the PY approximation. Santos et a/.0 
developed a rational function approximation to go beyond the PY approximation and de- 
rived improved expressions for the radial distribution functions and structure factors of SHS 
mixtures. 

Unfortunately, the PY analytical solution for mixtures requires the determination of 



a set of unknown density- dependent parameters A^, related, through p(p + l)/2 coupled 
quadratic equations, to other parameters which appear in the potentials as monotonically 
increasing functions of the temperature T and whose inverses measure the degree of adhesion 
( "stickiness" ) of interacting spheres. In most cases the coefficients Ay for given r%j can only 
be found numerically, and this feature limits the applicability of the SHS-PY model to small 
p values. As a matter of fact the number of actual applications to polydisperse fluids is very 
limited. We are aware of a study by Robertus et aZ@ on small angle x-ray scattering from 
microemulsions, with polydispersity represented by p = 9 components, a work by Penders 
and Vriji on turbidity of silica particles, and an investigation by Duits et on small 
angle neutron scattering from sterically stabilized silica particles dispersed in benzene. To 
simplify the numerical determination of the set {A^}, all these papers treat the special case 
of a single stickiness parameter, r, independent of particle size. 

In general, the SHS-PY solution for mixtures does not have the dyadic structure which 
allows the analytic inversion of Q (k) required to get closed expressions for structure factors 
of polydisperse systems. To recover the dyadicity and obtain an explicitly solvable model, 
Herrera and Blumil used the ad hoc assumption Ay = in a study on polydisperse CHS 
with sticky interactions under a MSA/PY closure. 

On the other hand, apart from Baxter's original definition^, there exists a second version 
of the SHS model, proposed in the one-component case by Brey et a/.0 and Mier-y-Teran et 
aZ.il. Here, the adhesive part of the potential is defined as the limit of a Yukawa tail when 
both amplitude and inverse range tend to infinity, with their ratio remaining constant. The 
analytic solution was obtained within the MSA closured, and Ginoza and Yasutomiil dis- 
cussed its relationship to Baxter's PY solution. Recently, Tutschka and KahlEl investigated 
the multi-component version of this second SHS model and presented MSA expressions of 
structure functions for both the discrete (finite p) and the continuous (p — > oo) polydisperse 
case. These authors pointed out that the dyadic structure of the SHS-MSA solution can 
be ensured a priori by imposing from the outset a Berthelot-type rulei on the coupling 
(stickiness) parameters 7^ of the Yukawa tail, i.e. 7y = (7^ 7^ ) 1 / 2 . 



Within the same framework of Yukawa-MSA models, the present paper has a threefold 
aim: i) we shall show that a new choice of Yukawa coupling parameters, Y i: slightly differ- 
ent from Tutschka and Kahl's ones0 ; can produce an even simpler solvable model of SHS, 
with a remarkable simplification of all analytical results; ii) we shall obtain closed analytical 
expressions for scattering intensity and other "global" structure factors of SHS by extend- 
ing the formalism successfully employed for HS and CHS; iii) we shall present numerical 
applications not only in the case of equal stickiness for all particles, but also when a simple 
size-dependence of the stickiness parameters Yi is assumed. 

The interplay among stickiness attraction, size polydispersity, and hard core repulsion 
gives rise to a rather complex behaviour. Nonetheless we shall present a simple unified 
description of these results, hinging on the introduction of a generalized Boyle temperature. 

The paper is organized as follows. In the next section the basic theory on structure 
factors, integral equations and dyadic matrices will be briefly recalled. The SHS model, its 
MSA solution under the assumption of factorizable coefficients, expressions for scattering 
intensity and other "global" structure factors will be given in Sec. III. Numerical results 
are included in Sec. IV, while the last section is devoted to a summary and some conclusive 
remarks. 



II. BASIC FORMALISM 

A. Structure factors and scattering intensity 

All scattering functions of multicomponent or polydisperse fluids with spherically sym- 
metric interparticle potentials can be expressed in terms of "partial" structure factors Sij(k), 
such as the Ashcroft-Langreth ones! 

S i:j (k) = Sij + Hij (k) = 5 i:j + {piPjf^hij (k) . (1) 

Here, k is the magnitude of the scattering vector, 8^ the Kronecker delta, pi the number 
density of species i, hij (k) the three-dimensional Fourier transform of the total correlation 
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function, hij (r) = (r) — 1, with gij (r) being the radial distribution function between two 
particles of species i and j at distance r. 

The knowledge of the Sij (k) allows to calculate the scattering intensity as well as some 
"global" structure factors. The coherent scattering intensity I(k) for a p-component fluid is 
given by&@ 

p 

R(k) =I{k)/V = p J2 (W) V2 F i (*) F I (*) Sij (k) , (2) 

where V is the volume, p = J2 m P m ^ ne total number density, while Xi = pi/p and Fi(k) 
denote the molar fraction and form factor of species i, respectively (the asterisk means 
complex conjugation). The measurable average structure factor is then defined from the 
Rayleigh ratio R(k) a ME 

S M (k) = R(k) j [ pP{k)\ . (3) 

with P(k) = Ylm=i Xm \ F m {k)\ 2 ■ As a third useful quantity, we consider the Bathia-Tornton 

number-number structure facto^, which is related to number density fluctuations: 

p 

<Snn (AO = ( x * x i) 1/2 • ( 4 ) 

The definition of other global structure factors may be found in Ref. 3. Clearly, R(k), Su (k) 
and S*nn (k) involve a unique kind of weighted sum, i.e., 

J2 «X (k) w* (k) Sq (k) , (5) 
with Wi (k) being equal to p l J 2 Fi (k) , [x i /P(k)] 1 ^ 2 Fi (k) , and x^ 2 , respectively. 



B. Integral equations in Baxter form 

Integral equations of the liquid state theory represent a powerful theoretical tool to get 
the (r) required to calculate the partial structure factors (k). The Ornstein-Zernike 
(OZ) integral equations relate the hij (r) functions to the simpler direct correlation functions 
Cij (r) . For fluids with spherically symmetric interactions, these equations arei@ 
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p , 

hij (r) = dj (r) + 22 Pm / dr' c im (r') /i mj - (|r - r'|) (6) 

m=l 

and can be solved only when coupled with an approximate second relationship (a "closure" ) 
among Qj (r), (r) and interparticle potential tHj (r)ii. 

By Fourier transformation, the OZ convolution equations become, in fc-space, 

[I + H (*)] [I - C (k)] = I, (7) 

where (A;) = {p i pj) l ^ 2 'cij (k) and Cjj (k) is the Fourier transform of Cy (r) . If S (/c) denotes 
the symmetric matrix with elements Sij (k), then we get 

S(k) = l + B.(k) = [l-C(k)]- 1 , (8) 

with I being the unit matrix of order p. The (k) can therefore be expressed in terms not 
only of hij (k), but also of (k) . However, in this paper we shall use a third representation 
of (k) based upon the Baxter factor correlation functions g^(r)I. By means of a Wiener- 
Hopf factorization of I — C (k) , Baxter transformed the OZ equations for HS fluids into an 
equivalent, but easier to solve, formi. Later on these equations were extended by Hiroike to 
any spherically symmetric potentials, without using the Wiener-Hopf factorization!. Baxter 
factorization reads 

I - C(k) = Q T {—k) Q (A;) , (9) 

where Q (k) has the form 

/+oo 
dr e ikr Q (r) , (10) 
-oo 

with Qij(r) = 2ix (piPj) 1 ^ 2 qij(r) ( Q T is the transpose of Q ). Note that Qij (—k) = Q*j (k) . 
For fluids of particles with spherically symmetric interactions including HS repulsions (i.e., 
Uij (r) = +oo when r < = (cr, + Oj-)/2, with <jj = hard sphere diameter of species i), the 
Baxter equations in r-space are 
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r Cij (\r\) = -Qijir) + 2vr J2 m Pm jT dt Qmi (*) (r + t) 



(11) 



rhij (|r|) = -%'(r) + 2vr £ m p m dt g im (t) (r - t) 

xh m j(\r-t\), 

where r > = (<jj — &j)/2 and the prime denotes differentiation with respect to r. 

Using Eqs. fl8|) and @, we get S(k) = Q 1 (k) Q 1 (— fc) and the partial structure 
factors can be written as 



S l3 {k) = £ (k) Q£ (-k) = Q7l (k) (k) 



On defining 



S U*0 = ]£>*(*)&J (*) 



i=i 



all the "global" structure functions can then be expressed as 

p p 
^ Wi (k) w* (k) Sij (k) = s m (fc)s^(/c). 

i,j=l m=l 



(12) 



(13) 



(14) 



C. Dyadic matrices and Vrij's summation 



The main problem of these analytical calculations hinges on the inversion of the matrix 
Q (k) — I — Q (k), which usually becomes a formidable task with increasing the number p 
of components. In a particular case, however, the inverse Q -1 (k) can be easily found for 
any size of the original matrix. This occurs when Qij(k) is a dyadic (or Jacobi) matrix, i.e. 
when it has the peculiar mathematical structure 



Q*i = Sa + i2 c ^ )b f (i,i = i,...,p) 



(the dependence on k was omitted for simplicity). We recall that a matrix TJ- 



(15) 



de- 



formed by the direct product of two vectors is often referred to as a dyad, ab, while a linear 
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combination of dyads ^2 A^a^b^ is called a dyadi^. Moreover, we shall refer to a sum 
of n dyads as a n-dyadic. 

We caution the reader that, since Qij(k) = 5{j — 2ir {pipj) %j (k) , where %j (k) is the 
unidimensional Fourier transform of (r) , Eq. fll5|) actually requires that qij (k) is a n- 
dyadic matrix (of order p), but in the following we shall use the same terminology for Qij(k) 
as well. The dyadicity is actually present in some solvable models of fluid mixtures: qij (k) is 
2-dyadic in the PY solution for neutral HstS, and 3-dyadic in the MSA solution for CHS0. 

The dyadic matrices have some special properties, which have already been partially 
discussed in Ref. 8. Here we recall the main points along with new additional features. Let 
us associate to the matrix of order p of Eq. (|15|), Q = I + Y^=i a^b^, a matrix Dq of 
order n with elements 

d a p = 8 aP + aW • (a, [3 = 1, . . . , n) , (16) 

where the dot denotes the usual scalar product of vectors, i.e. a^^ ■ b^= Y7 m =x a m°m- A 
first property is that any n-dyadic matrix Q of order p has always rank n, irrespective of p. 



Moreover, its determinant 



Q 



which is of order p, turns out to be equal to the determinant 
Dq = |Dq| of order n (with n <C p in multi-component fluids). A second property yields 
the explicit form of the elements of the inverse matrix Q _1 , as 

n n 
Q a=l /3=1 

where |DQ| a/3 is the cofactor of the (a,/3)th element in Dq and the double sum contains n 2 
determinants of order n. Clearly, this expression could also be rewritten as a sum of only 
n determinants, i.e., Q^ 1 = 8^ - YT a =i a i D j I 'Aj 5 where = Y?p=i°f IDqI"^ is the 
determinant obtained from Dq by replacing the a-th row with , ...,fi^ • Alternatively, 
one could write Qr. 1 = £J =1 Df ] bf ] /D Q , with d\ 0) = £™ =1 af ] |Dq|^ being obtained 
from Dq by replacing the /3-th column with af~\ (Cramer rule). 

All above expressions reflect the remarkable fact that Q _1 is a n-dyadic matrix as well, 
and it is indeed this property enabling a successful outcome of Vrij's summation for "global" 
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structure functions. Our starting point is a reformulation of Eq. ([L7]) in terms of a determi- 
nant of order n + 1 : 



1 



5ij 



(2) 



(I ; 



(n) 



6« 

1 + a« • b« 

a (2) . b (D 

aW • b« 



a (D . b (2) 

1 + a< 2 > • b< 2 > 



6} 



(n) 
3 

a« . b (n) 
a (2) . b (n) 



a («) . b (2) . . . 1 + a (») . b W 



(18) 



where Dq is included as a submatrix. From Eqs. ([T3|) and (|IBD, one then gets 



1 

A 1 



w ■ a^ du d\2 



b {2) ■■■ b (n) 



'22 



^2n 



(19) 



w • a (n) d nl d n2 ■■ ■ d nn 

To perform the sum over m required in Eq. (|H|), we expand this determinant along the first 
row to get 



3 m 1Um 



(20) 



where C a = T a /Dq and T a (a = 0, 1, . . . , n) is the cofactor of the element (1, a + l)th in 

the determinant of Eq. (|l9|). Clearly, T = -Dq. Using Eqs. fll4|) and (^0|), the searched final 
result is 

p n n n 

WjSij = w • w* + 2 ^ Re [w* • b^Cj b(Q) ' h((3) * C " C h ( 21 ) 

i,j=l a=l a=l (3=1 

where Re [...] denotes the real part of a complex number. VrijEi first performed a similar 
computation and derived a closed expression for the scattering intensity of HS mixtures. 
Our expression generalizes Vrij's one as well as that of Ref. 8 for the scattering intensity of 
CHS. It is simpler and more compact and can be used to calculate any "global" structure 
function, for mixtures with any number p of components. It involves a sum of only (n + l) 2 
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terms, and depends on p through some averages represented by scalar products of vectors. 
Only the number of terms contained in these averages increases with increasing p, and hence 
the application to polydisperse mixtures is straightforward. 



III. THE MODEL 

A. Sticky hard spheres as a limit of Yukawa particles 

A fluid of SHS can be derived from particles interacting via HS plus Yukawa (HSY) 
attractive potentials, i.e. 

{+00, < r < an 

(22) 

in the limit p — > +00 with /3Aij = pKij and independent of p. Here, f3 = (fcgT) -1 , where 
ks is the Boltzmann constant, and all Aij are positive, with Aji = Aij and Kji = Kij, as 
required by the symmetry condition Uji(r) = Uij(r). This approach is convenient, since the 
Baxter equations for HSY mixtures have been solved analytically!!, for any finite p, within 
the MSA closure, which adds to the exact hard core condition, hij (r) = — 1 for r < a^, the 
approximate relationship Cy (r) = —f3uij{r) for r > a^. The solution is 



0, r < Lij 



\<n{r - a^) 2 + (bi + ai«7y)(r - <r y ) + Qj [e^-^ - l] 
+D i3 e-^ r -^\ <r< a i3 



(23) 



Dije-^ r -^\ r > a tj 



where the coefficients are determined by a complicate set of equations^. From these, how- 
ever, it can be shown that, as p — > +00, then Cy — > —D%j and Dy — > (3 Aij/ p = K^. The 
MSA solution for SHS is therefore 

\ai(r - a^) 2 + (6; + Oi<7y)(r - <r y ) + K ij} Uj < r < a tj 
qij{r) = <( (24) 

0, elsewhere 
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V V 

^ = \Y. P**?' Ci = 1 , A = 1 - 6- (26) 
i=i j=i 

To calculate Q'^k), we need the unidimensional Fourier transform ^ (A;), i.e. 

AX, ll,,^., J_„ 3 aw) 



^ 1 + ^_ 12Ci J_,|^ (27 ) 

where X m = ka m /2, jo(x) = sinx/x and ji(x) = (sinx — xcosx) jx 1 are spherical Bessel 
functions, and i - when it is not a subscript - is the imaginary unit. 

B. Factorizable coefficients 



In general q^ (k) , as defined in Eq. (|2T[), does not have the required dyadic structure, 



due to the presence of in the last term. To overcome this difficulty, Tutschka and Kahl0 
proposed the following Ansatz: 

Kij = lij a ij, wit h lij = {la Ijj ) 1/2 (Berthelot-rule), (28) 

which yields = 7* ; 7// 2 (of + 2<Ji<jj + er|)/4. Consequently, the last term of Eq. ([27D 
splits into three independent contributions, and g^- (/c) turns out to be 5-dyadic, in spite of 
the fact that in the HS limit (no adhesion) it is only 2-dyadic. 

We first note that a great simplification occurs with the factorization 

= YiYj (29) 

(all Y m > 0). In this case, the last term of Eq. (E7|) generates only one contribution, and 
%j (k) becomes simply 3-dyadic and, in a particular case to be discussed later on, even 2- 
dyadic. The defined by Eq. ( |2"9"D satisfy the Berthelot-rule, i.e. = (KaKjj) 1 ^ 2 . Note 
that the stickiness parameters 7 mrn are dimensionless@, while the Y m are lengths. 
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Factorizable adhesive parameters have already been considered by Yasutomi and 
GinozaH and Herrera et a/.0 in studies on adhesive-HSY fluids, although no expressions 
for structure functions were given. Since in those papers = GGiGj, the relationship 
with our coefficients is simply given by Y m = G l l 2 G m . 

Using Eq. fl2Tj|), we get Q = Yi with = (tt/6) Y7j=i Pj a j^j (dimensionally analogous 
to £2), and therefore 

— - 1/2 

Qij(k) = Sij - 2n (pipj) %j (k) 

= ^^{(i + ^-i%-^r t )&*& (30) 



1/2 
3 ' 



which has the required dyadic structure Q%j = + YlZ=i a i ■ We emphasize that the 
decomposition into and bj is not unique. Our choice allows an easy comparison with 
the corresponding results for polydisperse Hsi and CHSi. After defining; 



a m{k) = ttt ' 



m) 

m 



2A m X m 

= ^J (X n ), (3 m (k) = + - i0 a m (k), (31) 

6£\k) = -1-no m Y m j Q (X m ), 8 m (k) = $£>(fc) - I2g a m (k), 

Qij(k) may be rewritten as 

Q..(k) = 5 tJ + p} /2 e iX > [aj(k) + cnPjik) + YiSjik)] p)'\ (32) 

and the corresponding decomposition is 



(33) 





= p!'V* 


»r 


1/2 
= Pj <*j, 




= Pl'^-a, , 




l/2n 
= Pj Pj, 


«f 




6f 


l/2r 

= Pj S 3- 



1 /2 

If we are interested in the scattering intensity, then Eq. (19) with w m = pm F m yields 
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F 

- 1 m 


Cm 


/3m 


&m 


1/2 
Pm 


{F} 


i + M 


m 


{S} 


Dq 


{aF} 


{era} 


1 + {a(3} 


{a5} 




{YF} 


{Ya} 


{Yf3} 


1 + {Y5} 



(34) 



where Dq coincides with the cof actor of the (l,l)th element in the 4x4 determinant, and 

p v 
{fg} = p m e iX -f m g m = P (e iX fg) , (f) = ^ x m f m . (35) 



m=l 



m=l 



Here and in the following, angular brackets (• • •) denote compositional averages over the 
distribution of particles (this notation differs from that of Ref. 8, which also contains 
misprints corrected in Refs. 35,36). 

With the chosen decomposition of Qij(k) the stickiness contributions are confined in the 
last row and the last column of the determinant of Eq. ([34]) . When adhesion is turned off, all 
these elements vanish apart from their diagonal term, and the 4x4 determinant essentially 
reduces to the 3x3 HS one. For numerical computation, it is convenient, following Ref. 
8, to simplify s m by using elementary transformations which do not alter the value of the 
determinant. If we add to the third column the second one multiplied by — (3^/A — ik/2) 
and to the fourth column the second one multiplied by 12^, then s m becomes 





F 

- 1 m 




1/2 
Pm 


{F} 


1 + {«} 


D Q 


{aF} 






{YF} 


{Ya} 



0, 



(0) 



(36) 



{/?(°)} -3&/A + ifc/2 
{Y(3^} 

where Dq has been changed accordingly. Expanding the 4x4 determinant along the first 
line and inserting it into Eq. (0), the final result for the scattering intensity of SHS is 



{c7(°)} + 12£ : 
{a8®} 
1 + {Y5^} 



R{k)/p = (F 2 ) + (a 2 ) |Q| 2 + </?(°) 2 > \C 2 \ 2 + (S^ 2 ) \C 3 \ 2 
+2 Re [(Fa) C x + (F(3^) C 2 + <F5 (0) ) C 3 
+ (a/3<°>) C X C* 2 + (aS®) dC* + (/W°>) C 2 C 3 *] 



(37) 
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where form factors have been assumed to be real quantities, as is indeed the case for spherical 
homogeneous scattering cores. On the r.h.s. all k arguments have been omitted for simplicity, 
and the C v {k) have already been defined with reference to Eq. (|20|). 

The expression for the average structure factor Sm (k) is then obtained after division of 
the Rayleigh ratio R(k) by pP(k) = p (F 2 (k)) . Moreover, the Bathia-Thornton number- 
number structure factor Snn {k) can easily be derived by setting all F m = 1 everywhere into 
the expression of R(k) / p, with the result 

Snn (k) = 1 + p «« 2 > |d| 2 + </^ 2 ) \C 2 \ 2 + <^ 2 ) |C 3 | 2 

+2 Re [(a) C x + </3 (0) > C 2 + (<5 (0) ) C 3 (38) 
+ (a/3®) C x Cl + (aS^) C X C* + (^S^) C 2 Q] } , 

where the C u are the analogues of the C u appearing in R{k). 



IV. RESULTS FOR POLYDISPERSE FLUIDS 
A. Size distribution 

For SHS fluids containing only one chemical species, size polydispersity simply means the 
presence of a multiplicity of possible diameters. In a "discrete" representation of polydisper- 
sity the number p of different diameters is very large but finite, and Xi is the fraction Ni/N of 
particles having diameter <jj. On the other hand, a theoretical representation with infinitely 
many components (p — ► oo) and "continuously" distributed diameters is also possible and 
often used. 

Although all formulas of previous Sections refer to a finite number p of components, 
the polydisperse continuous limit of such expressions can immediately be inferred by the 
replacement rules x a — > dx = f(a)da and "$2 x a ... — ■> J daf(a)..., where f(a)da is the 
fraction dN/N of particles with diameter in the interval (a, a + da). As molar fraction 
density function f(a) we choose the Scholz distribotionBi 
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/((x) = r4y ( a u) exp (" a M) (a>1) ' (39) 

where T is the gamma function!!, (a) the average diameter, a = l/s 2 , and s a = 
[(°~ 2 ) — ( cr ) 2 ] 1 ^ / ( a ) measures the degree of size polydispersity. In the monodisperse limit, 
s a = 0, the distribution becomes a Dirac delta function centered at (a). The Schulz function 
allows an easy analytic evaluation of some averages J daf(a)..., such as the moments (cr m ), 
which obey a simple relation for m > 1, i.e. 

m— 1 

(O = [1 + (m - 1)4] (a) = J] M, , (40) 

with Mj = 1 + jsl . 

In most cases, however, analytical integration is hardly feasible, and numerical integra- 
tion brings back to discrete expressions with large p, of order 10 2 — 10 3 . In practice, the 
"discrete" representation of polydispersity is the most convenient for numerical purposes, 
and all formulas of the previous Sections can be employed by assuming x a = f(a a )Aa, 
where Aa is the grid size of numerical integration. 

For fluids with Schulz-distributed diameters the packing fraction, r] = £ 3 = (n/6)p (a 3 ), 
can be written as i] = r] mono (1 + s 2 a ) (1 + 2s 2 ) , with ?7 mono = (vr/6)p (cx) 3 . 



B. Stickiness distribution 

On a dimensional basis, the parameters must be lengths. Moreover, = YjYj must 
be proportional to j3 — (ksT)^ 1 . If we assume, for simplicity, that stickiness polydispersity 
and size polydispersity are fully correlated, then the most natural choice for Yj is 

Yi = 7o <?i , (41) 
with the dimensionless proportionality factor 
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where Eq denotes an energy and Yu\ is the stickiness parameter of particles with diameter 
(a) . This implies that 

Kij = 7q <JiOj = OiOj = ^ OiOj , (43) 

where we have also introduced a reduced temperature T* = (k B T/e ) = I/Jq. 

The model of Eq. fl4"T| ) will be compared with the one of SHS polydisperse in size but 
not in stickiness (on the analogy of Refs. 20 and 21). In this simpler case all particles have 
the same Fj = = 7 (a) , and the degree of stickiness polydispersity sy, defined similarly 
to s a , vanishes. 

Both these models may be regarded as particular cases (for a = and a = 1) of a more 
general size-dependence given by 

^=^(m)" =7 °^ (44) 

with a > 0. We have examined this generalization for a = 2 and a = 3, but for the purposes 
of the present paper we restrict our analysis only to the cases a = and a = 1. 

The choice Y{ = 70 <Ji has very interesting properties. First, the corresponding distribu- 
tion of F-values, related to the size distribution j a as fy = dN/dY = f a da/dY, is a Schulz 
function as well, with (Y) = Y/ a \ and Sy = s a . A second more important fact is that only 
in this special case Qij(k), in general 3-dyadic for SHS, becomes simply 2-dyadic, i.e. 

Qa{k) = S tj + p] /2 e iX > {A-(*0a;(*0 + G a4°\k)}p) /2 , (45) 

with 

AW = 1 + (^ - i|) „ , 
G„ = 1 - 4 1o 2 A = 1 - [1 - , mono (1 + si) (1 + 2sl)\ (46) 



(now £2 = 70 ^2)- It is remarkable that this expression for Qij(k) differs from the HS one 
only for the presence of Gq (Gq = 1 for HS). Now the natural dyadic decomposition becomes 



„(1) _ 1/2 iX t A. , (1) _ 1/2 

a i -Pi e A, Oj- - Pj a„ 

„(2) _ l/2.iX, r _ l(2) _ _l/2fl(0) ' 
a i — Pi e <^0&i > Oj — Pj Pj , 

while the 4x4 determinant appearing in Eq. ( |34]) reduces to a 3 x 3 one, with a conse- 
quent simplification of the formulas for R(k)/p, Sm (k) and 5"nn (k) (all terms depending on 
subscript 3 vanish in Eqs. (^) and Q3"8P). 

C. Numerical results 

Because of its importance in the analysis of experimental scattering data, we have 
focused on the measurable average structure factor Sm (k). The scattering cores inside 
the particles have been assumed to be spherical and homogeneous, with form factors 
F m = V^ catt An m 3ji(X^ catt )/X^ catt , where X^ catt = A;< att /2, < att < a m is the diameter of 
a scattering core of species m, V^ catt = (ir/6) (o"m att ) 3 ^ s volume, and An m its scattering 
contrast with respect to the solvent. For mixtures with several components belonging to 
only one chemical species, as in the present paper, An m is the same for all particles. For 
simplicity, we have taken a^ att = a m . 

The polydisperse SHS model depends on the following parameters: the packing fraction 
rj, the strength 7 of the adhesive interaction, the average diameter (a) , and the two degrees 
of polydispersity s a and sy- In all numerical calculations we have adopted dimensionless 
variables, with lengths expressed in units of (a) . To understand the influence of each param- 
eter on Su(k), it is instructive to first recall the behaviour of a sequence of simpler systems, 
starting from monodisperse hard spheres and adding in the first two cases either surface 
attraction or size polydispersity. 

1. Monodisperse HS and SHS 

In pure fluids all particles are equal (s a = = sy), r\ = r] mono , and Sm (k) = S mono (k) 
with no form factor involved. Figures 1 and 2 depict the dependence of S mono (k) on the 
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parameters (rj, 70). Fig. 1 illustrates the evolution of S^ono {k) , as 77 increases from low values 
up to the freezing one, in the well known case of monodisperse HS of diameter a without 
stickiness (70 = 0). Here we have exploited the PY solutioni, which, for HS, coincides with 
the MSA one. In Fig. 2 the dependence of Smono (k) on 77 is displayed for monodisperse SHS, 
at two fixed 70 values, i.e. 0.5 and 0.7, corresponding to 7q = e / (fc^T) = 0.25 and 0.49, or 
to T* = 4 and 2. 04 (70 and T* have been defined in Eqs. (f42[) and (]43|), respectively). 



In all these cases (Figs. 1, 2a and 2b), as r] increases at fixed 70, the first peak height and 
amplitudes of all subsequent oscillations increase, but the behaviour near the origin depends 
on 70, as will be discussed in more detail shortly. 

On the other hand, the effect of increasing 70 (i.e. increasing the adhesive attraction or 
decreasing T) at fixed r\ can be seen by comparing, for instance, the solid curves (77 = 0.2) 
of Figs. 1 and 2. As 70 increases, the first peak and subsequent maxima are shifted to 
larger k values, and their amplitudes change as well. However, the most significant effect 
on Smono (k) occurs near the origin. Here, S mono (0) substantially increases and becomes the 
global maximum at large 70. This behaviour can be understood from the explicit expression 
of S mono (0), which reads 



Smono (0) = Qmono (0) = " — — -3 (48) 

J [1 + 277 - 12 7 2 77 (1 - 77)] 
= Kt/K' 1 * = pk B TK T . 

Since S mono (0) is related to the isothermal osmotic compressibility K T and to the density 
fluctuations^, its drastic increase signals the approach to a gas-liquid phase transition. The 
critical point can be obtained from the spinodal line, defined by S~l no (0) = 0, and the 
critical parameters turn out to be: rj c = (\/3 — l) /2 ~ 0.37 and 7q c = (y/3 + 2) /6 ~ 0.620 
(corresponding to 7o c — 0.79, or to the reduced critical temperature T* ~ 1.61). 

The combined influence of rj and 70 can be observed going back to Fig. 2. On defining the 
Boyle temperature Tg as the one where the attractive and repulsive forces balance each other 
in such a way that the second virial coefficient B2 vanishes^, we note that the temperatures 
corresponding to 70 = 0.5 and 0.7 lie, respectively, above and below Tg (but in both cases 
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above T*). In fact, for this monodisperse model it is easy to see, from the low-density 
expansion of S mono (0) = 1 + ph (0) ~ 1 - 2B 2 p + C(p 2 ), that B 2 = AV m {\ - 3 7o 2 ) = 
4Vks(l-3/T*) with Vks = (^/6)ff 3 , and therefore T* = 3 (7 2 B = 1/3 or 7oB ^ 0.58). 
Fig. 2 suggests the existence of two different "regimes" for Smono (k) above and below the 
Boyle temperature, respectively. When T* > Tg or, equivalently, 7 < 7ob ("weak-attraction 
regime" , as in Fig. 2a) the fluid behaves like pure HS without stickiness. Here repulsive forces 
are dominant, B 2 > 0, and compressibility and density fluctuations, along with the whole 
^mono (k) near the origin, decrease with increasing rj. When T* < T* < Tg or, equivalently, 
7ob < 7o < 7oc ("strong-attraction regime"), one finds B 2 < 0, while the balance between 
attractive and repulsive forces becomes more complex. In this case Smono (k) near the origin 
has a non-monotonic dependence on rj, as in Fig. 2b. Here, compressibility and density 
fluctuations first increase with rj, in agreement with the low-density expansion of S m0 no (0). 
Then, an inversion occurs at 770 = (6 — 2T*) / (6 + T*) (~ 0.24 when T* = 2.04) and afterwards 
Smono (0) decreases. In other words, below Tg attractive forces seem to be dominant at low 
packing fraction, whereas repulsion again prevails at higher rj. 



2. Polydisperse HS without stickiness 

Fig. 3 refers to polydisperse HS without surface adhesion (70 = 0). Size polydispersities 
s (j = 0.1,0.3 have been employed here and in the following, since values in this range are 
rather common in experimental data from colloidal fluids. The two Schulz distributions 
have been discretized with a grid size Act/ (a) = 0.02, and truncated where f(a)Aa ~ 10~ 8 , 
i.e. at cr cut / (a) = 1.68 and 3.48, respectively. Since each diameter characterizes a different 
component, these discrete polydisperse mixtures involve p = 85 and 175 components. Note 
that these numbers of components are much larger than those used with the SHS-PY model 
of Ref. 20. 

The effect of size polydispersity is considerable@, as appears from a comparison among 
Figs. 1, 3a and 3b: with increasing s a at fixed 77, Sm (k) slightly increases in the low-fc region, 
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its first peak is reduced and shifted to smaller k values, and all subsequent oscillations are 
progressively dumped, as a result of destructive interference among the several length scales 
involved. 

3. SHS polydisperse in size but not in stickiness 

At this point we study SHS fluids polydisperse in size but monodisperse in stickiness, 
with all particles having Yj = = 7 (a) (s CT 7^ 0, a = =>- s Y = 0). This choice will be 
referred to as Model I and has been prompted by the SHS-PY investigations of Refs. 20-22, 
where a single, size-independent, stickiness parameter was considered. 

Figures 4 and 5 illustrate what happens when a surface adhesion (with 70 = 0.5,0.7) is 
added to size polydispersity. Comparison with Fig. 3 (70 = 0) shows that the attractive 
interaction, in the presence of size polydispersity, produces a further lowering of oscillation 
amplitudes in the first peak region and beyond. When 7 = 0.7 and s a = 0.3 (Fig. 5b) all 
curves exhibit an almost complete flattening in the same range. 

Near the origin (for k (a) < 5), for both considered cases with 70 = 0.5 (weak- attraction) , 
only a small increase in Sm (k) is found with respect to polydisperse HS without stickiness 
(Fig. 3), the relative ordering of all curves is unchanged and also coincides with that of the 
corresponding monodisperse SHS (Fig. 2a). On the other hand, when 70 = 0.7 and s a = 0.1 
(strong-attraction and low size polydispersity, Fig. 5a) the behaviour of Sm (k) close to the 
origin strongly differs from that of polydisperse HS without stickiness and is similar to the 
monodisperse SHS case of Fig. 2b. Surface adhesion produces large Sm (0) values, which are, 
however, smaller than the corresponding monodisperse ones. This means that, when 7 7^ 0, 
size polydispersity reduces Sm (k) even near the origin (whereas, when 70 = 0, increasing 
s a at fixed rj determines an increase of Sm(0)). This effect of size polydispersity, in the 
presence of attraction, is amplified when 70 = 0.7 and s a = 0.3 (strong-attraction and high 
size polydispersity, Fig. 5b). Now one observes an interesting return to a "HS-like ordering" 
of the curves in the low-fc region, as in the case 70 = 0.5. This behaviour is peculiar of 
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Model I and will be absent in Model II to be presented in the next subsection. 

4- SHS poly disperse both in size and in stickiness 

Next we consider the case of stickiness correlated to the size, according to the linear law 
Yi — 70 (Tj (a — 1, sy — s a ^ 0). This will be referred to as Model II. 

The results for 70 = 0.5 , shown in Fig. 6, are qualitatively similar to those of Model I 
(Fig. 4). When sy = s a = 0.1 the quantitative differences are very small. However, when 
Sy — s a — 0.3 the Sm(0) values lie more clearly above those of Fig. 4b. 

For 7 = 0.7 (Fig. 7) the behaviour of Sm (k) in the first peak region and beyond is 
essentially unchanged with respect to Model I, but near the origin differences are larger and 
significant. Here, when sy = s CT = 0.1 the Sm (k) curves are similar to those of Fig. 5a, 
with larger Sm (0) values (very close to the corresponding monodisperse ones of Fig. 2b), 
but as sy = s a = 0.3 there is a qualitative as well as quantitative difference with respect 
to Model I (Fig. 5b). Indeed in the low-fc region the Sm {k) curves of Fig. 7b exhibit 
the same relative ordering present in the previous case with lower polydispersity (Fig. 7a) 
as well as in the corresponding fully monodisperse fluid (Fig. 2b). This persistence in a 
"strong-attraction regime" even at high size polydispersity constitutes the main difference 
between Model I and II. Such a feature can be probably related to the fact that the stickiness 
distribution of Model II is skewed towards large Yi values completely absent in Model I, 
and this asymmetry implies, on average, stronger attractive forces. 

Unfortunately, the behaviour of Sm(0) in polydisperse models does not admit any simple 
thermodynamical interpretation. For mixtures, in fact, the average structure factor Sm {k) 
depends on the form factors and Sm (0) is no longer the normalized compressibility!!. Nev- 
ertheless, we have been able to account for the aforesaid difference of "regimes" between 
Model I and II when 70 = 0.7 and s a = 0.3 in terms of a single parameter, which generalizes 
the Boyle temperature of the monodisperse SHS case. 
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5. Generalized Boyle temperature 



The Boyle temperature of these polydisperse models can be found by deriving their 
second virial coefficient B 2 from the low-density expansion of Snn (0) = 1 — 2B 2 p + 0(p 2 ). 
Likewise, to interpret the behaviour of Sm (k) previously discussed, we start from the low- 
density expansion of Sm (0) . A straightforward calculation, employing the dyadic formalism 
of Sections II and III and not reported here, yields 

§y X (0) = - p{x i x j ) 1 '' 2 [| (<rj + 3a,a 2 ) - ItsY&p^ + 0(p 2 ), (49) 



and therefore, from Eq. fll2|), 

Sij(0) = Sij - pixiXj) 1 ' 2 [| (a, + a,) 3 - 2tt (FjY^fTj + ^Ifo)] + 0{p 2 ). (50) 

Inserting this result into Sm (0) = J2i,j( x i x jY /2 [ F i (°) *7 (°) / p (°)] % (°) > and usin S the 
above mentioned expression for F m (k), we obtain Sm (0) = 1 — 2B 2F p + ^(p 2 ); with 

«V = i^5) K" 6 > <- 3 > + 3 {"') ~ 12 (° 3 Y) <^'>] • (51) 

which is the analogue of the second virial coefficient, including all form factors. The sign of 
B 2F , and therefore the behaviour of Sm (0) at low density (as well as the overall "regime" 
of <Sm(&))> depends on T*, which is hidden in the Y terms. On defining a generalized 
Boyle temperature Tg F as the one where B 2 p vanishes, and employing our assumption 
Yi = Jotf/ (o-)"' 1 , we find 

T* = U{a3+a) (g4+Q) (52) 

B ' F [(<7 6 )(<7 3 ) + 3(cr 5 >(c7 4 >] (a} 2ia - 1} ' 

Exploiting Eq. fl40|), it follows that 

Tg F = 12/ [M 4 (M s + 3M 3 )] for Model I, 

(53) 

T£ F = 12M 3 / (M B + 3M 3 ) for Model II, 
where Mj has been defined with reference to Eq. (flop. The role of Tg F for Sm (k) is the 
same as that of Tg for S^ono (k) of monodisperse SHS: above Tg F the behaviour is "HS-like" , 
whereas a "strong- attraction regime" is found when T* < T* < Tg F . 
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Eqs. ( |53|) imply that, in both models, Tg F depends on the degree of size polydispersity 
s /j . However, whereas Tg F (s CT ) of Model I is a rapidly decreasing function, in Model II 
it exhibits only a very slow decrease asymptotically approaching 18/7 ~ 2.57. Such a 
difference explains the behaviours displayed in Fig. 5 and 7, which refer to T* = 2. 04 (i.e. 
7o = 0.7). In Model I, when s a = 0.1 one has T* < T£ F = 2.79, whereas when s a = 0.3 
one finds T* > Tg F = 1.68. On the other hand, in both cases the temperature of Model II 
is below the generalized Boyle one, the values of Tg F now being 2.99 and 2.90 for s a = 0.1 
and 0.3, respectively. 

Finally, it is instructive to compare F with the true Boyle temperature Tg of these 
polydisperse models. The second virial coefficient, obtained from the low-density expansion 
of S NN (0) , turns out to be B 2 = (vr/6) [(a 3 ) + 3 (a 2 ) (a) - 12 (Y) (aY)] , and one obtains 
Tp^ = 12/ [Mi (M a + 3)] and 12/ (M 2 + 3) for Model I and II, respectively. Note that for these 
polydisperse fluids it is always Tg < (T^) mono = 3, T^ F < Tg for Model I, and T£ F > for 
Model II, while in the limit of monodisperse fluids Tg F — > (Tg) mono . 

V. SUMMARY AND CONCLUSIONS 

In this paper we have presented a new analytically solvable model for multi-component 
SHS fluids within the MSA closure, using a hard-core- Yukawa potential with factorizable 
coupling parameters (in the appropriate infinite amplitude and zero range limit). The model 
is simpler than previous ones available in the literature, since qij(k) is in general 3-dyadic 
(Tutschka and Kahl's mo del§ was 5-dyadic), with a consequent great simplification of all 
analytical formulas. 

We have stressed the importance of the "dyadic structure" of qij(k) and recalled the 
properties of matrices with dyadic elements. Such a feature allows the analytic inversion 
of Qij(k) required to get the partial structure factors Sij(k). Through Vrij's summation, 
expressions have then been obtained for global structure functions, such as R(k), Sm (k) and 
Snn (k). These closed analytical formulas, just as their counterparts for polydisperse HsEi 
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and CHS0, allow to "bypass" the computation of the individual p(p+ l)/2 partial structure 
factors, which may be a rather difficult task for polydisperse systems with large number 
of components. Because of their simplicity, our expressions may therefore represent a very 
useful tool to fit experimental scattering data of real colloidal fluids. 

While the presented 3-dyadic expressions hold true for any choice of stickiness param- 
eters Yi, two particular versions of the model have been analyzed numerically. The first 
one assumes size polydispersity, but a single stickiness parameter for all particles (Model I), 
while the second one proposes stickiness parameters dependent on the diameters according 
to a linear law (Model II). Model I is similar to the SHS-PY models for polydisperse col- 
loids known in the literature&~§l, while Model II is the simplest choice for size-dependent 
stickiness parameters. 

The combined influence of hard core repulsion, adhesive attraction and polydispersity can 
generate a variety of behaviours at the level of measurable average structure factor Sm (k). 
We have recognized the existence of two different "regimes" for Sm (k) both in monodisperse 
and in polydisperse SHS fluids. Above a temperature Tg F , which in the monodisperse case 
coincides with the Boyle temperature, we have identified a "weak-attraction" behaviour, 
resembling the HS one. In the range below Tg F but still above the critical temperature, a 
"strong-attraction" regime sets in, and we have described its features in detail. It is found 
that Tg F is a decreasing function of the degree of size polydispersity s a . It is also worth 
noting that the behaviour of our SHS-MSA models in the "strong-attraction regime" is in 
qualitative agreement with that of the SHS-PY model displayed in Fig.s 3 and 4 of Ref. 20, 
where the existence of two different "regimes" for Sm (k) was, however, not recognized. 

All our numerical results show that size polydispersity strongly affects the behaviour of 
Sm {k) in the first peak region and beyond, where the influence of stickiness polydispersity 
is less significant. Models I and II are nearly equivalent in this interval of k- values, whereas 
they may substantially differ near the origin. The present study shows that the use of 
a single stickiness parameter, instead of more realistic size-dependent ones, may lead to 
marked differences in the small angle scattering region at sufficiently high 70, i.e. when 
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attraction is strong or temperature is low. 

In the small k region the adhesive forces and the specific relationship between sticki- 
ness and size parameters have far reaching consequences. Although very little is known 
experimentally about the correlation between stickiness and size, it is reasonable to expect 
that larger particles attract each other more strongly. The linear dependence Yj = 70 <Ji 
represents the simplest non-trivial choice, but other possibilities could also be taken into 
account@. Duits et a/.@ found that in some cases the SHS-PY model with a single sticki- 
ness parameter, independent of particle size, was unable to fit their experimental scattering 
data, and these authors already emphasized the possible role of stickiness polydispersity as 
a cause for the observed deviations. Similar discrepancies between experimental and model 
Sm (k) values could be a crucial test for the soundness of our Model II with respect to Model 
I, as well as of any other choice for the stickiness-size functional relationship. 

It would be also instructive to compare our SHS-MSA model with other recent theoretical 
approaches to polydisperse colloidal fluids. As an example, we mention the "optimized 
random phase approximation" joined with orthogonal polynomial expansions, proposed by 
Lado and coworkers^. 

However, the most important advantage of the present model lies in its simplicity. In 
particular, version II has special formal properties, since - only in this case - the expression 
of qij(k) becomes 2-dyadic. Therefore version II can indeed be reckoned as the simplest 
solvable model for polydisperse SHS and it could be a good candidate to tackle the issue of 
thermodynamics and phase stability of these fluids from a fully analytical point of view. We 
expect that compact expressions for pressure, chemical potentials, partial structure factors 
at k = 0, as well as other quantities required to investigate - for instance - sedimentation^, 
vapor-liquid equilibrium and demixing in the presence of polydispersity, can easily be ob- 
tained. We hope to accomplish this task in a forthcoming paper. 
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FIGURES 

FIG. 1. PY (= MSA) structure factor S mono of monodisperse hard spheres without stickiness, 
plotted as a function of the dimensionless variable ka at various packing fractions ij. Note that 
this and all following figures have the same scale. 

FIG. 2. MSA structure factor 5 mon o of monodisperse sticky hard spheres, plotted as a function 
of ka at various packing fractions 77, for two values of adhesive strength parameter 70 (here, and 
in all following figures, the curves of Part b) are shifted upwards by 2.5 units). 

FIG. 3. PY (= MSA) structure factor Sm of polydisperse hard spheres without stickiness, 
plotted as a function of k (a) at various packing fractions rj, for two different degrees of size 
polydispersity s a . 

FIG. 4. MSA structure factor Sm of sticky hard spheres polydisperse in size but monodisperse 
in stickiness (Model I), plotted as a function of k (a) at various packing fractions rj, for 70 = 0.5 
and two different degrees of size polydispersity s a . 

FIG. 5. As Fig. 4 but for 7o = 0.7. 

FIG. 6. MSA structure factor Sm of sticky hard spheres polydisperse in size and stickiness 
(Model II), plotted as a function of k (a) at various packing fractions rj, for 70 = 0.5 and two 
different degrees of size polydispersity s a . 

FIG. 7. As Fig. 6 but for 7o = 0.7. 
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